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Abstract. Fractional Hamiltonian Monodromy is a generalization of the no- 
tion of Hamiltonian Monodromy, recently introduced by N. N. Nekhoroshev, D. 
A. Sadovskii and B. I. Zhilinskn for energy-momentum maps whose image has 
a particular type of non-isolated singularities. In this paper, we analyze the 
notion of Fractional Hamiltonian Monodromy in terms of the Gauss-Manin 
Monodromy of a Riemann surface constructed from the energy-momentum 
map and associated to a loop in complex space which bypasses the line of sin- 
gularities. We also prove some propositions on Fractional Hamiltonian Mon- 
odromy for 1 : — n and m : —n resonant systems. 

1. Introduction 

We consider an integrable system on a four dimensional symplectic manifold 
defined by an energy-momentum map. For a proper map, the Liouville-Arnold 
theorem allows to foliate the phase space by tori or a disjoint union of tori over the 
regular values of the image of the map. Hamiltonian monodromy is the monodromy 
of this fibration |Dui80|, ICB97] . The word Hamiltonian is added to distinguish 
this monodromy from the Gauss-Manin monodromy of Riemann surfaces which is 
also used in this paper. A non trivial monodromy can be expected if the set of 
regular values of the image of the energy-momentum map is not simply connected. 
Hamiltonian monodromy has profound implications both in classical and quantum 
mechanics |Ngoc99 since it is the simplest topological obstruction to the existence 
of global action-angle variables |Dui80| and thus of global good quantum numbers 
|Ngoc99| . The phenomenon of Hamiltonian monodromy has been exhibited in a 
large variety of physical systems both in classical and quantum mechanics AKE04 , 
ISCOOl IKR031 ICWT991 E7M IW.IDnSl IGCSZ04I [EJSM] , 

The presence of non-trivial monodromy in energy-momentum maps with isolated 
singularities of focus-focus type is now well-established. The non-trivial monodromy 
in the spherical pendulum is due to this singularity F)ui80| ICB97] . Recently, the 
definition of Hamiltonian monodromy has been extended to characterize not only 
isolated singularities but also some types of non- isolated singularities, leading to 
the concept of Fractional Hamiltonian Monodromy |NSZ02| INSZ061 IEfs04j [ECSu7] . 
More precisely, one considers an energy-momentum map with a 1-dimensional set 
C of weak critical values defined by the property that each point of this set lifts 
to a particular type of singular torus, a curled torus, i.e. for the simplest case two 
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cylinders glued together along a line whose extremities are identified after a half- 
twist. The standard Hamiltonian monodromy describes the possible non-triviality 
of a 2-torus bundle over a loop in the set of regular values of the image of the 
energy-momentum map. The monodromy matrix is an automorphism of the first 
homology group Hi of the torus with integer coefficients. Fractional monodromy 
can appear if the set of admissible paths is enlarged to include loops which cross 
the singular line C. For such paths, the singularity of C being sufficiently weak, 
it can be shown that the monodromy action can still be defined but only on a 
subgroup of Hi. The formal extension of this action to the whole group leads to 
monodromy matrices with fractional coefficients and to the denomination fractional 
monodromy. One of the main motivations for the introduction of this new concept 
is given by the quantum manifestation of monodromy in the discrete joint spectrum 
of the energy- momentum map [Ngoc99, NSZ06J. This spectrum can be represented 
as a lattice of points in ]R 2 . The focus-focus singularity can be detected by a point 
defect of this lattice which prevents it to be a regular lattice isomorphic to 1? . 
In the same way, fractional monodromy can be interpreted as a line defect of the 
lattice and appears therefore as a natural generalization of standard hamiltonian 
monodromy. 

The presence of fractional monodromy has been shown in a system of coupled 
oscillators in m : —n resonance with morn different from 1. Two constructions have 
been given based on geometric [NSZ06, Nek07 or analytic [ECS07J arguments to 
define rigorously the crossing of C. The geometric construction consists in following 
a basis of cycles of a regular torus through the crossing of C. Not all the cycles 
can cross continuously the singularity, only those corresponding to a subgroup of 
Hi can. When C is crossed, one allows cycles to break up and reconnect, the 
orientation of the cycles being preserved. The second construction uses, as in the 
original paper of Duistermaat [Dui80, CB97 , the period lattice of the torus [Arn89 
which is however not defined on the singular line C. The period lattice is defined 
through two functions and T at each point of the regular values of the image 
of the energy-momentum map. Some regularizations of these functions can be 
made in order to cross continuously the line of singularities |ECS07| . Note that the 
preceding geometric point of view can be reconstructed from this analytic approach 
since the basis of cycles can be determined from the functions and T. 

The goal of this work is to study fractional monodromy by complexifying the 
phase space in order to bypass the line of singularities. Somewhat similar studies 
for the Lagrange top and the spherical pendulum have already been published 
|Aud021 IViv031 IBCOlj and have highlighted the relat ion between Hamiltonian and 
complex monodromy. A parallel can also be made with Bohr-Sommerfeld rules 
for semi-classical quantization. Such calculations can be undertaken in the C°° 
[dVP99| |NgocOO[ ldVN03j or in the analytic context |Vor83| IDDP971 IDP97] . The 
real approach needs regularization of the sub-principal term whereas the complex 
approach avoids such problems by avoiding the singularity. In this paper, we show 
that fractional hamiltonian monodromy is given by a Gauss-Manin monodromy 
AGZV88, IZ0IO6] of a Riemann surface constructed from the energy- momentum 
map. The construction can be made in the reduced phase space [CB97] which is 
well suited to the introduction of Riemann surfaces. The Gauss-Manin connection 
is defined for a complex semi-circle around the line C of singularities and can be 
calculated by applying the Picard-Lefschetz theory. We also introduce the complex 
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extension of the functions and T which are viewed as integrals of rational forms 
over a cycle of the Riemann surface. In other words, the regularizations of and 
T in the real approach are replaced by a complex continuation of these functions. 
The variations of and T along the bypass in the complex domain are deduced 
from the Gauss-Manin monodromy. The functions and T allow us to go back to 
the real approach and to define a real monodromy. We show that this monodromy 
corresponds to fractional hamiltonian monodromy asymptotically in the limit where 
the radius of the complex semi-circle goes to 0. Using this construction, we recover 
the results of the real appro ach obta ined in Refs. [NSZ02llNSZ06UEfs04UECS07] for 
the l:-2 resonance and Rcf. Nek07 for m : —n resonance. Moreover, for 1 : — n and 
m : — n resonant systems, we give new proofs of these results from the real and the 
complex approaches. A geometric point of view of the Gauss-Manin monodromy 
can be given by inspecting the motion of the branching points of the Riemann 
surface along the semi-circle around C. For 1 : —n resonant systems, the Riemann 
surface has locally n complex branching points in a neighborhood of C lying on a 
circle around the origin with an angle of 2-7T jn between each other. Along a complex 
semi-circle around C, the n points turn by an angle of 2ir/n and exchange their 
positions. The variation of near C is then calculated as a residue of a given 
1-form. This characterizes the line of singularities C and fractional monodromy. 

The organization of this article is as follows : We first consider the real approach 
and we determine the monodromy matrices for l:-2, 1 : — n and m : — n resonances 
in Sec. O We next show in Sec. [3] how fractional monodromy can be defined in 
the complex approach and we detail its geometric interpretation. We recover the 
different results obtained in the real approach. Concluding remarks and perspec- 
tives are given in Sec. 2J Appendices O and [B] present a schematic representation 
of the geometric construction of fractional monodromy for l:-2 and 1 : — n resonant 
systems in the real approach and the semi-classical point of view for 1 : — n and 
m : —n resonances. The material of these appendices complements the existing lit- 
erature on these two points. Appendix [Cl finally deals with the reduction procedure 
in the complex approach. 

2. The real approach 

The goal of this section is to furnish a short overview of fractional hamiltonian 
monodromy in the real approach. Starting from the example introduced in Ref. 
Efs04. ECS07 for the l:-2 resonance, we extend it to 1 : — n and m : — n resonances. 
The core of the results presented in this section are already contained in [NSZ02, 
INSZ061 IEfs04| IECS07] . We describe it in some detail since we need the results 
for the complex approach. Furthermore, we present analytic computations in this 
section and a geometric construction in appendix which are slightly different 
from the original ones and give a new view on fractional hamiltonian monodromy. 
The originality of the analytical computation for the l:-2 resonance relies in the 
introduction of a local description of the energy-momentum map near the origin 
of the bifurcation diagram which considerably simplifies the computation of the 
monodromy matrix. Since these arguments are generalizable to 1 : — n and m : 
— n resonances, they allow us to prove some propositions stated in NS Z061 IEfs041 
Nek07 . For the geometric description, we introduce the standard representation of 
a torus i.e. a rectangle whose some edges are identified. This geometric construction 
can be extended straightforwardly to 1 : — n resonant systems. 
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2.1. The l:-2 resonance. We consider the symplectic manifold M = T*R 2 with 
standard symplectic form lu = dqi A dpi + dq2 A dp2- We introduce the energy- 
momentum map F = (H, J) : M — > R 2 where the two functions J and H have zero 
Poisson brackets {J,H} = 0. Following Refs. [NSZ021 INSZ061 lEfsOl IECS07| . we 
choose a system corresponding to the 1 : — 2 resonance defined by 

m / H = V2[(q 2 1 -p 2 1 )p2 + 2q lPl q 2 }+2e(q 2 1 +p 2 1 ){q 2 2 +p 2 2 ) 

where e is a non-zero real number. We denote by 1Z the image of F and by 7Z reg 
the regular values of TZ. We recall that a point M £ 1Z is regular if the 1-forms dH 
and dJ are linearly independent at all points of F (h,j). 

The flow of J defines an 5 1 -action on the phase space but this action is not prin- 
cipal since the isotropy groups of the points {p\ —Q,q\ — 0,p2,<?2} are isomorphic 
to Z/2Z. The reduced phase space J~ 1 (j)/S 1 can be constructed by using the alge- 
bra of invariant polynomials with values in R which is generated by |Efs04|, ICB97] 



(2) 



J(p,q) = l/2 (q!+pi)~(q 2 +p 2 2 ) 
7Ti(p,q) = 1/2 (qj+pj) + {ql+p 2 2 ) 
7T2 (P,q) = \/2[(?i -pf) q_2 - 2qipxp 2 ] 
7r 3 (p, q) = V2[(qf - pf) p 2 + 2gipiQ 2 ] 



with the constraint \J\ < tti. The reduced phase spaces Pj = J 1 (j)/S 1 are defined 
in the space R 3 = (tti, ^2,^3) by the equations : 

(3) ^+^. 2 = ( 7ri _ j) ^ 1+J ) 2 ; 

and correspond to non-compact surfaces with a conical singularity for j < 0. Having 
introduced the invariant polynomials, some comments can be made on the choice 
of H. Since {H,J} = 0, it can be shown that if H is polynomial in (pi,qt) then 
it can be written as a polynomial function in J, 7Ti, 7r 2 and 7r 3 . We also notice 
that the reduced phase space for the l:-2 resonance being non-compact, not every 
choice of H leads to a proper map for the energy-momentum map F — (H,J). 
This explains why a term of degree four has been added to H [see Eq. (JIJ)]. The 
image of the energy-momentum map defined by Eq. {T]) has the particularity to 
present a line of singularities C, each point of this line except the origin lifts to a 
singular curled torus (the origin lifts to a pinched-curled torus). The topology of 
these singular tori can be determined by the intersection of the reduced phase space 
Pj with the level set Hj = h (see Sec. 12. 4p . Hj is the reduced Hamiltonian, i.e., 
a map from Pj to R that sends a point of Pj to H(tti,tt2,tt3, j). Moreover, since 
the two energy- momentum maps F — (H, J) and F' = (H — /(J), J), where / is a 
polynomial function, define up to a diffeomorphism the same fibration of the phase 
space, we can consider an example such that H = for the line of singularities C. 
Figure Q] displays the bifurcation diagram of F (see |Efs04[ IECS07] for details on 
this construction). Singular points (i.e. not regular) are represented by solid lines 
in Fig. [TJ C is a line of weak singularities. The word weak means that for each 
point of C, there exist points of the corresponding pre-image such that the rank of 
F is 1. We recall that for a regular point, the rank of F is always equal to 2 and 
that the rank of F is for one point of F^ 1 (Q, 0). 
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Figure 1. Image 1Z of the energy-momentum map F (in grey). 
The singular line C is represented by the horizontal solid line. The 
small full dot indicates the position of the image of the pinchcd- 
curled torus. A loop T crossing transversally C is depicted in 
dashed lines. 



2.2. Fractional Hamiltonian Monodromy. We begin this section by recalling 
some basic facts about integer and fractional monodromy. The word hamiltonian 
in hamiltonian monodromy will be omitted when confusion is unlikely to occur. 

Standard or integer monodromy is defined for a loop T along regular values of 1Z. 
Fractional monodromy is obtained by extending the possible loops, allowing them 
to cross some particular type of weak singular lines such as C. In the standard 
case, we consider a point (h, j) of T which lifts to a torus T 2 (h, j). We fix a basis of 
the homology group H\(T 2 (h, j),Z). There exists a natural connection of the torus 
bundle which allows to transport this basis along T [CB97, Bat91 . The monodromy 
matrix, which is an automorphism of Hi , is the holonomy of this connection. This 
construction can be generalized to fractional monodromy but only a subgroup of 
H 1 (T 2 (h,j),Z) can be transported continuously across the line C |NSZ061lECS07| . 
These remarks can be understood by the following construction. Let (h,j) G H r eg- 
Denoting by tpj and <ph the flows associated to the Hamiltonians J and H, the 
period lattice of F at a point (h, j) is the set 

(4) {(ti,t2)eRV/° = 

for all z 6 A basis for this period lattice, which is isomorphic to Z 2 , 

is given by the vectors V\ = (27r,0) and = (— 0,T) where is the rotation 
angle and T the first return time of the flow tp H defined as follows Dui80, CB97J. 
The Hamiltonian J generates an S^-action on j). We denote by 9 an angle 

conjugated to the action J. Following cpj which is parameterized by 9, one goes 
back to the starting point when 9 increases by 2-7T, which gives v\. If we consider 
now ipH from a point of an orbit of the flow ipj, one sees that the first intersection 
of these two flows takes places at time T. The two points of intersection of the two 
flows define two angles 9/ and 9i and the twist = 9 f — 9i which is determined 
with respect to the direction of ipj. Note that a different choice of the angle 9 leads 
to a different basis for the period lattice. The corresponding rotation angles 
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(G) 



Figure 2. Schematic representation on the torus F 1 (h, j) of the 
cycles fa and fa depicted respectively in dashed and solid lines. 

differ by a multiple of 2ir. The monodromy matrix associated to a loop lying in the 
regular values of the image of the energy-momentum map is related to the behavior 
of the functions O and T along this loop. For the standard monodromy, after a 
counterclockwise loop around an isolated critical value (focus-focus singularity) it 
can be shown that the rotation angle is increased by 2tt whereas the first return 
time is unchanged. As a consequence, v\ is transformed into v\ and v 2 into — v\ +i>2, 
and thus the monodromy matrix M written in the local basis (vi,V2) is equal to 

( 5) m=(_\ ; 

The two functions 8 and T allow to define a basis of cycles for the homology group 
Hi(F^ 1 (h,j),'Z) and thus to recover a more geometric point of view. A basis 
([Pi], [fa]) of j),Z) is given by the cycles associated respectively to the 

flows of the vector fields 

Xi = 2nXj 

X 2 = -Q(h,j)Xj + T(h,j)X H ■ 

The flows ip Xi (t £ [0, 1]) and ip X2 (t 6 [0, 1]) generate respectively the closed cycles 
fa and fa. These two cycles are schematically represented in Fig. [2] for O = — ir. 
In the basis [fa]), the monodromy matrix is given by the same matrix as Eq. 

§5§ obtained for the period lattice. 

The situation is slightly more complicated for fractional monodromy. The ana- 
lytic construction of fractional monodromy follows the same steps as for the stan- 
dard case }ECS07j . Returning to the example of Eq. ([I]) and considering the loop 
r of Fig. [1] the question which naturally arises is the definition of the crossing of 
the line C since O has a discontinuity of size ir on this line and T diverges. The 
idea of the method proposed in Ref. ECS07 consists in prolonging by continuity 
the function (the continuous function is called 0) and in rescaling the time to 
obtain a finite first return time denoted r. Note that this rescaling does not modify 
the definition of the cycles but only the time to cover them. This leaves therefore 
the monodromy matrix unchanged. O is defined by Q(h,j) = Q(h,j) for points 
(h, j) of r before the crossing of the line C and by Q(h,j) = Q(h,j) + tt for points 
after |ECS07j . A basis of the period lattice is given by the vectors v\ = (2-7T, 0) and 
«2 = (— ©, r) and the cycles (i\ and fa are now associated to the vector fields 

X 1 = 2ttX.j 

x 2 = -e(h,j)Xj + T(h,j)x H ■ 



(7) 



The construction of these cycles has however to be carefully examined. ip f x (t £ 
[0, 1]) generates a cycle fa for all points of T. Before the crossing, <p^ (t £ [0, 1]) 
generates a closed cycle fa but after the crossing, ip^ (t £ [0,1]) generates only 
half of a cycle. To get a complete cycle, we thus have to take t £ [0, 2]. This means 
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geometrically that only a cycle /3 2 covered twice can be transported continuously 
across the line C. Thus, not all cycles in the homology group can be transported 
along r. Only a subgroup corresponding to the cycles that are run twice by X 2 can 
be transported. In terms of the period lattice, the crossing of C is only possible for 
the sublattice generated by v\ and 2v2- We can then define the monodromy matrix 
for a counterclockwise loop crossing transversally C once. The monodromy matrix 
M reads in the basis (v\, 2v 2 ) or in the basis [2/?2]) 



(8) M 



1 

-1 1 



Extending formally the definition of M to the whole homology group or the whole 
period lattice, we obtain in the basis (^1,^2) or in the basis ([/3i], {P2]) 



(9) M 



1 



We finally note that this generalized monodromy is still topological in the sense 
that it depends only on the homotopy equivalence class of the loop considered. We 
also point out that these cycles prolonged continuously to the curled torus allow to 
recover the geometric construction of fractional monodromy described in appendix 

El 

2.3. Local computation of the monodromy matrix. The determination of the 
monodromy matrix is based on the behavior of the functions and r on the line of 
singularities C. More precisely, the monodromy matrix can be constructed uniquely 
from the size of the discontinuity of and from the fact that r is continuous. In 
Ref. |ECS07| , the computation was done by using global expressions for and t in 
terms of elliptic integrals. It is clear that such a global calculation can be expected 
to be done explicitly only for simple energy-momentum maps. 

We propose a computation of the fractional monodromy matrix based on lo- 
cal arguments for the energy- momentum map F, i.e., a Puiseux expansion in h 
and j around (h = 0,j = 0). This expansion is not trivial because a particular 
dissymmetry in h and j has to be preserved. 

Lemma 1. In the variables (tti, 7r 2 , 7r3, J) and for a point (h,j) € TZ r eg, the func- 
tions and t are given by the following expressions : 

(10) &{h,j)=h 



G?7Tl 



% U +7Tl)v /< 3( 7r l) 



(•/// 



d 



(11) r(h,j) = l f 1 J ifi=dn 1 , 

1 j*z v < 9( 7r i) 

where Q(tti) is a polynomial given by 

(12) Q(7n) = (7n-j)(7n+j) a -[fc-e(7r?-i 2 )] a ■ 

■Ki and TTi are the two largest real roots of Q with ttJ~ < tt^ . 

Proof see |ECS07j and Sec. (|2.5|) for a more general construction. We recall that 
t is the first return time of the flow of the rescaled vector field tUIh.D 
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Proposition 1. The monodromy matrix for a counterclockwise oriented loop T 
crossing once the line C transver sally at a point different from the origin ( see Fig. 
QP is given by 



(13) M 



1 

-1/2 1 



Proof The monodromy matrix is given by the behavior of and r in the neigh- 
borhood of the line C. Taking j < fixed and finite, the limits lim^Qi Q(h, j) and 
lim/ w0 ± r{h,j) have been calculated in [ECS07 . Elliptic integrals and asymptotic 
expansions of these integrals were used. 

We propose a simpler computation by considering the asymptotic limit j — > 0. 
For that purpose, we analyze the roots of the polynomial Q as h and j go to zero. 
Since a qualitative change of the functions and r is expected when the polynomial 
Q has a multiple complex root, we determine the complex discriminant locus of Q 
near the origin h = j = 0. This point will be made clearer with the introduction 
of Riemann surfaces in Sec. [3] but here it gives the way in which the two limits 
h — > and j — > should be taken. The discriminant locus in the real approach 
has already been calculated since it corresponds to the line of singularities of the 
bifurcation diagram, i.e., to the points where the 1-forms dH and dJ are linearly 
dependent (see Fig. [1]). We introduce the variable x = j +7Ti and h, j and x are 
taken complex. Three roots of Q vanish for h = j = 0. 

Let us assume that x, j and h go to zero. Constructing the Newton polyhedron 
associated to Q |Kir93j . we obtain the principal part Qn of Q which can be written 



2 



(14) Q N = x 3 - 2jx 

We notice that this principal part is symmetric with respect to h which is not 
the case for the polynomial Q. Simple algebra leads to the following asymptotic 
discriminant locus 



(15) 



h = ±^P(-§) 
h = 




denoted A. This complex locus is displayed for j € K in Fig. [51 Equation (TT4")) also 
shows that the weights 2, 2 and 3 can respectively be attributed to x, j and h. In 
other words, if we introduce the small parameter r, the principal parts of x, j and 
h can be written 

(16) 

The monodromy can be computed along a small loop in 1Z around the origin. This 
loop can be parameterized by (h,j) £ S 1 which allows to derive a local version 
of the computation of fractional monodromy. We thus let h — > while keeping j 
fixed and finite. Note that it is equivalent to consider the limits h,j—+0 with the 
condition h = o(j 3 / 2 ), since on A we have asymptotically h — 0(j 3 ^ 2 ). In order to 
determine the leading terms of the roots of Q in this case, we construct the Newton 
polygon of Qn. The approximate solutions fulfill 

(17) - 2jx 2 - h 2 = . 
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Since the sum of the three roots of Q that go to zero when h,j—>0 is 2j, and the 
sum of the four roots is , one deduces that the principal parts of the roots are 
given in the original variables by 



(18) 



2j 

-h 
7? 



These expressions can be compared with the ones given in Ref. }ECS07j where 
j < is fixed and h — > : 



2j 

-h 



(19) 



-2j+2ej 
h 



On _|_ 1 _|_ Mj. 

47 + p + J+2^J 



We now calculate lim^ ± j <0 Q(h, j) and lim^ j <g T(h,j). We consider first 
h > and j < 0. Q(h,j) can be written as 

rXi dx 

(20) ' ' 



h r 

*e Jx 3 



"y{x — x\)(x — X2){x — xz)(x — Xi) 



We determine only the principal term of the asymptotic expansion of 0. The 
symbol ~ represents the equivalence in the limit h — ► 0, j — ► and h = o(j 3 ^ 2 ). 
From Eqs. (fl8|) . we obtain 



(21) Q(h,j)^- [ 



h/V=2J xV(x - 2j)[x + h/y/-2j)[x - h/^-2j){x - l/e a ) 
We decompose the preceding integral into three integrals by introducing the terms 
k and fe' which go to such that J^} 2 j "C & <C — 2j <C A;' <C 1- fc and k' are chosen 
for instance as j a . The notation a <C & means that the ratio a/6 — > as a and 6 
go to 0. The three integrals are taken over the intervals [\h\/y/—2j, k], [k,k'] and 
[k' , 1/e 2 ]. It can be shown that the limit of the last two integrals is zero. The first 
integral reads 

(22) (->(/<../:•- - ' 

which can be rewritten as 



(23) &(h,j)~h 



h/V^J xy/(-2j)(x*-hy(-2j))(-l/e*) 
dx 



h/V =75 Xy /{-2j)(x*-hy(-2j)) 



The change of variables x = . „ . u leads to 

+oo 



du 



u^u 2 — 1 



(24) Q(h,j)r 

Using the fact that 

(25) / — dU = arctan[V" 2 - 1] 
J uyv? — 1 
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one finally arrives to 

(26) lim 0(fc,j) = E. 

h->0+ ,j— 0- ,h=o(j*/ a ) 2 

Similar calculations for h < give 

(27) lim Q( h ,j) = ^. 

We then deduce that the discontinuity of is equal to it. 

t can be calculated along the same lines. This time the first two terms go to 
zero and we only determine the last one 

1 f 1 ^ 2 xdx 

(28) r(h,j) 1 



2«£ Jk' y/{x - Xi)(x — X2){x — X S )(x — X4) 

Simple algebra leads to 

1 f 1 /^ dx 
(29) r(h,f 1 

and using the fact that 



2*£ J k , ^x(x ~ 1/e 2 ) ' 

(30) / , . d \, 2 =2\n{^+^x—W 2 

J \/x(x - 1 e 1 ) 



y/x(x - 1/e 2 ) 
one obtains that 

(31) lim r(h,j) = . 

t is therefore continuous on the line C. The monodromy matrix is finally deduced 
from the behavior of and r in the neighborhood of C. We follow for that purpose 
the construction of Ref. [ECS07] which is briefly recalled in Sec. 12.21 □ 

Remark 1. The preceding computation being local does not show the topological 
character of fractional monodromy i. e. its independence with respect to the homo- 
topically equivalent loops considered or more simply with respect to j. This point 
has been proved in the real approach in Ref. [ECS07] and will be proved in the 
complex approach in Sec. 

2.4. 1 : — n resonant system. We consider the energy-momentum map F — 
(H, J) where J is given by 

(32) J=\[(vl + ql)-n(ql+pl)], 

with n > 2. We assume that the bifurcation diagram of F is locally in the neighbor- 
hood of the origin given by Fig. [3] The singular locus corresponds to (h = 0, j < 0). 
Each point of this locus lifts to a singular torus i.e. a pinched-curled torus for the 
origin and an n— curled torus for the other points. A fc-curled torus is a singular 
torus for which one cycle is covered k-times while the others only once. We now 
restrict the discussion to a particular family of energy-momentum maps having lo- 
cally the bifurcation diagram of Fig. [3l As in Sec. 12.21 the first step consists in 
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h 



Figure 3. Local bifurcation diagram in the neighborhood of the 
origin (h = 0, j = 0) for the 1 : — n resonant system. The singular 
locus is represented by the large solid line. The small full dot 
indicates the position of the origin. The dashed line depicts a loop 
used to calculate the fractional monodromy matrix. 



determining the invariant polynomials associated to the momentum J. We have 
[Efs04] 



(33) 



J(p, q) = l/2[m (q\ +p\)-n (q 2 2 + pf)] 
TTi (p, q) = l /2[m (g f + p\) + n (g| + pf)} 
7T 2 (p, q) = Vn m m n ^[( qi + ipi) n (q 2 + if> 2 ) r 



?T3(p, q) = Vn m m n 5i[(qi + ipi) n (q 2 + ip2) m } 
for m > 1 and n > 1. The reduced phase space Pj is defined by 

(34) A+^l = ^i+j) n ^i-]) m , 

with the condition tt\ > 

Definition 1. FFe consider the set T of energy-momentum maps F — (J, _ff ) which 
can be written as 



(35) F 



J 

H = tt 3 +R(it 1 ,J) 



where R is a polynomial so that F is a proper map. ~K\, tt^ and J are given by Eqs. 
(Mfy form = 1. 



Note that we do not search to determine or characterize the set T. Only some 
properties of the elements of T will be sufficient to compute the monodromy matrix. 
Simple examples can be exhibited to show that T is not empty. For instance for 
the resonance l:-3, we can choose 



(36) F 



J 

H = 7r 3 -(7r 1 - J)(tti + J) 4 



We are interested in the local behavior of R near the origin or in other words under 
which conditions on R, the image of the corresponding energy-momentum map is 
given by Fig. [3l 



12 



D. sugny\p. mardesic 2 ,m. pelletier 2 : a. jebrane 2 ,h. r. jauslin 1 



Lemma 2. The energy-momentum map given by Eqs. LS5\) has locally the bifurca- 
tion diagram of Fig. [3] in a neighborhood of the origin if the polynomial R is of the 
form 

(37) Rfa , J) = (7T1 + Jf (7T1 - J)™' fl(7Tl , J) , 

where n' and m! are positive integers such that n' > § and n' + m' > . R 
is a polynomial in tt\ and J such that the two smallest real positive roots of the 
polynomial Q which are larger than j for (h, j) G TZ reg are simple roots . Q is the 
polynomial defined by 

(38) Q = (n 1 +j) n (n 1 -j)-[h-R(TT 1 ,j)} 2 . 
The two roots are denoted tt^ and tt^ with ir^ < tt^ . 



Proof The proof is based on the nature of the intersection of the reduced phase 
space Pj with the level sets {Hj = h} of equations h = tt3+R(tti, j) as h and j vary. 
Fig. 2] displays these intersections for three different values of h, j < being fixed. 
The case considered in this figure is the l:-3 resonance and the energy- momentum 
map of Eqs. ([55)) . Fig. ^represents the generic topology of the level sets {Hj = h] 
which we are going to characterize. We recall that each point of the reduced phase 
space lifts in the original phase space to a circle except for the point of coordinates 
(ni = — j, 7T2 = 0,7T3 = 0) which lifts either to a circle covered n times, if j < 0, 
or to a point, if j = 0. One deduces from the bifurcation diagram of F that the 
intersection of the level set {Hj = h} with Pj contains a circle passing through the 
singular point S — (iri = —j, ir 2 = 0, 7r 3 = 0), for j < and h = 0. We set y = n 3 , 
x = it i + j and x 1 = tt\ — j . The local behavior near the point S of the two surfaces 
is given by 

(39) y = ±x%x'^ , 
for Pj and by 

(40) y = -R(7T 1 ,j), 

for {Hj — 0}. It is then straightforward to show that the local behavior expected 
is obtained if R(jri,j) — x n x' m R(TTi,j) with the conditions n' > Q and n' + 
m' > ^ji, m' > 0. The first and second inequalities result respectively from the 
conditions for j < and j = 0.D 

Remark 2. The inequalities of lemma&are strict to ensure that a multiplication by 
a constant factor of the term R does not modify the local behavior of the image of the 
energy-momentum map. This point has not been assumed for the l:-2 resonance, 
the parameter e is thus chosen sufficiently small in this case. Here, this hypothesis 
simplifies the computation of the monodromy matrix as can be seen in the proofs of 
lemma^ and proposition^ 



Proposition 2. The monodromy matrix of a l:-n resonant system of the form 
i35\) is given for a loop turning counterclockwise around the origin ( see Fig. and 
crossing once transversally the line C by 



(41) M 



1 

-1/71 1 
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Figure 4. Intersections of the reduced phase space Pj (in large 
solid lines) with the level sets Hj = h for h > 0, h = and h < 
in the plane 7T2 = 0. The energy-momentum map corresponding 
to this diagram is given by Eqs. (|36|) . 

Remark 3. Propositions [H and\3\ (see Sec. \2.5\) were formulated as conjectures 
in Refs. NSZ06, Efs04 . They are based on the analysis of the quantum joint 
spectrum of the energy-momentum maps. They have been recently proved in Ref. 
|Nek07j from a geometrical construction. Note also that our starting point here and 
in Sec. \2.5\ is more general than in Sec. \2.2\ in the sense that only a local structure 
of the bifurcation diagram of the energy-momentum map F = (H, J) is assumed. 

As was done for the l:-2 resonance, we have to determine the expressions of the 
functions and t in terms of the invariant polynomials, the monodromy matrix 
being given by the behavior of these two functions in the neighborhood of C. 

Lemma 3. The functions and r are given for a point (h, j) € IZreg by the 
following integrals 



Oo is a function over 1Z which is continuous on the line of singularities C and 
which therefore gives a trivial contribution to the monodromy matrix, t is the first 
return time of the rescaled vector field X — Xjj- 

Proof By definition, we have for a point (h, j) € IZreg that [ECS07] 



where 9 is the angle variable conjugate to J . It can be expressed in terms of the 
variables (<?i,Pi) by 



(42) 




(43) 




(44) 



6 = arg(pi +iqi) ■ 
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This corresponds to a particular choice of the angle 9. Other choices lead to the 
same monodromy matrix. This point will be detailed in Sec. 12.51 for the m : — n 
resonance. Differentiating Eq. (|44[) and using Hamilton's equations, one arrives to 

(45) 9 = dv /\ 9qi , 
which simplifies into 

a 1 , , n M OR . 

(46) 9 = — ( n h-nR+—pi + — qi ). 

J + 7Ti dpi dqi 

The integral of Eq. (|43[) can be rewritten as an integral in the reduced phase space 
Pj [ECS07| : 



(47) Q(h,j) = 2 



7Tl 



Using the particular form of the polynomial R (see lemma ^ , it can be shown that 
the last three terms of Eq. (|4"6"1) gi ye a continuous contribution to the function 
denoted 9o- Since tti = 2mi2 [Efs04 , we finally obtain that 

(48) &(h,j)= r j^^= + Q (h,j) . 

J«l J + TTl VQ 



The term tt2 has been replaced in Eq. (|48|) by combining Eqs. (|34|) and (|3G|) . 
The determination of r is straightforward if we remark that 

(49) r(h,j) = £ ds = ^dt = J\p\ + qlT^dt , 

where s and t are respectively the rescaled and the original time. The rest of the 
proof consists, as we did for O, in rewriting the integral of Eq. (j49|) in the reduced 
phase space Pj and leads to 

(so ^^r^,,,. 



□ 

The last technical point to be discussed is the behavior of the roots of Q as h go 
to zero. 

Lemma 4. The complex discriminant locus A of Q near the origin (h = 0,j = 0) 
is given by 



(51) A = 



h = 

^ h = ± \Z _n ' i ^" +1 7n+l)('«+ 1 ) 

The polynomial Q as a function of x and in the limit h — > 0, j < fixed has n roots 
Xk whose leading term reads 

(52) X k = a e 2<*fc/* 

with k S {0, 1, 2, n — 1}. TTie other roots of Q have a non- zero finite limit. 
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Proof We first determine the complex discriminant locus A of Q near the point 
(h = 0, j = 0). Constructing the Newton polyhedron of Q and taking into account 
only the terms of lower degrees, the principal part Qn of Q can be written 

(53) Q N (x) = x n+1 - 2jx n - h 2 . 

A straightforward calculation then leads to A. 

In the limit h — ► 0, j < fixed, we construct the Newton polygon associated 
to Q. The roots of the principal part of Qn satisfy —2jx n = h 2 which allows to 
deduce the n roots Xk (k S {0, 1, • ■ • ,n — 1})-D 

We have now all the tools ready to prove proposition [21 

Proof We first consider the case h > and j < 0. We recall some of the properties 
of the roots of the polynomial Q viewed as a function of x which will be used in 
the calculation. We assume that Q has N roots with N > n, denoted X{ (i € 
{0, 1,---N— 1}). x n ,x n +i, ■ ' ' iXn-i are the roots of Q of order 1 in the limit 
h — * 0, j < fixed. x n is the smallest real positive root of Q with a non-zero limit. 
Since the polynomial R is defined up to a multiplicative constant, we can write Q 
without loss of generality as follows Q(x) = YVi^o 1 ( x ~~ x i)- Examination of the 
coefficients of Q leads to the following relations 

f n:=r i ^=(-i) jv+1 ^ 2 

(54) mro^-i)-- 1 ^ ■ 
I nt; 1 ^=2j(-i) jv - n+i 

We determine only an equivalent of the function 8 — Bo- We obtain 

dx 

(55) e(h,j)-e (h,j) ■ ' 



i 2/ " l—ri=N 



oTS x 



TJ i=N-l, 

rL=o ( x 



wnere (-2jy/™ ^ s ^ ne smallest positive real root of Q. We introduce the function k 

such that (J_2j)i/" ^ k 1 when h — > 0. We then proceed as in the proof for the 
resonance l:-2 by decomposing the integral of Eq. ([55]) into two integrals. Only the 
first integral from rr^jjT/n to k has a limit different from zero. We then have 



(56) e(h,j)-G (h,j)~h 



The change of variables x = j^jy/^u leads to the following expression for — 6o 

[ + °° du 

(57) e(h,j)-Q (h,j) 



dx 



^4-2j)(x»-^)xl 



Vu n - 1 



Using the fact that 



du 2 



(58) / — . = - arctan[Vu n - 1] , 
J u^Ju n - 1 n 

one finally arrives to 

(59) Hp e(hj)-e (h,j) = - . 
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Figure 5. Same as Fig. |3]but for the m:-n resonant system. 



It can also be shown that 

(60) hm e(h,j)-e (h,j) = — . 

h— i-o-,j<o n 

Following the same arguments, we can calculate the limit of r. We decompose r 
into two integrals but we only examine the last one denoted r%. The calculation of 
the other integral can be done along the same lines. An equivalent of ti is given by 

1 x n ~ x dx 



(61) T 2 (hJ) 
which can be rewritten as 

1 f Xn dr 

(62) 75(ft,j)~r 



y/x- n + 2 (x N - n - 2j) 
which has the same finite and non-zero limit as h — > ± with j < fixed. □ 

2.5. Generalization to m : — n resonances. For the m : — n resonant case, the 
momentum J reads 

(63) J=\[™(p\+<ll)-n{ q l+pl)] , 

where m > 2 and n > 2 are relatively prime integers. The bifurcation diagram 
corresponding locally near the origin to an m : — n resonant system is displayed in 
Fig [5] The equation of the singular locus is h = 0. Each point of this line lifts for 
j < to an n-curled torus whereas the points of C with j > lift to an m-curled 
torus. The origin corresponds to a pinched-curled torus. 

Proposition 3. The monodromy matrix of a m : —n resonant system of the form 
i35\) is given for m and n relatively prime and for a counterclockwise loop around 
the origin (see Fig. by 



(64) M = 



1 
^ 1 

mn 
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The proof of proposition [3] follows the same lines as for the resonance 1 : —n. As 
in Sec. 12. A\ we consider the family of energy-momentum maps given by Eqs. (|35j) 
where tti, 1T3 and J are given by Eqs. (|33[) with m > 1. Some lemmas are required 
before the final proof. 

Lemma 5. The energy-momentum map F of Eq. H35\) has locally the bifurcation 
diagram of Fig. in a neighborhood of the origin if R is of the form 

(65) i?(vr 1 , J) = (tt! + J)"' (ttx - J) rn ' R(tti , J) , 

where n 1 > j, m! > ^ • R is a polynomial such that the two smallest real positive 
roots of Q which are larger than j are simple roots for (h,j) € TZreg- Q is the 
following polynomial 

(66) Q = ^i+3) n ^i-j) m -{h-R^ l ,j)) 2 . 
The two roots are denoted and . 

Proof The proof is similar to the proof of lemma |2JD 

Lemma 6. The functions and r are defined for a point (h,j) G !Z r eg by the 
following integrals 



^j) = ^/5 (J+ ' l) ""^ 1 ~ jr " 1 ^i 



(67) 

I ;-!/, ,'l — 

0o is a function over 1Z which is continuous on the line of singularities C . u and 

v are two integers such that mu — nv — 1 . r is the first return time of the rescaled 

vector field X = , ■> , -^^-i, ; , -i~, m -i Xh defined such that t has a non-zero finite 
(Pi+n) (P2+12) 

limit on C. 

Remark 4. One can obtain the expression given by Eq. | f^ff| ) corresponding to the 
case 1 : — n as a special case of Eq. with u = 1 and v = 0, which is a particular 
solution of mu — nv = 1 for m = 1. 

Proof The determination of O is based on the dependence of the angle 8 as a 
function of the coordinates (pi, qi,P2, 12)- To clarify this question, we introduce the 
canonical conjugate coordinates (Ii,ipi) and (12,^2) which are defined as follows 



(68) 



qk = V^h sin ip k 
p k = cos ip k 

Note that the polar coordinates (Ik,<fik) are only defined if p\ + q? > |Arn89j 
which is not the case on the line C. By definition, the angles ipi and tf2 vary 
in an interval of length 2n. We look for a linear canonical transformation which 
transforms the two angles ipi and tf2 into 9 and ip where the angle is canonically 
conjugate to J. The angular dependence of H in the variables (Ji, <pi, I2, ^2) is 
given by the term ir^ and is equal to rupi + mip2- We then set 



(69) 



V = Utpi + Vif2 

V> = nipi + mip2 



where u and v are integers such that mu — nv = 1 which ensures that the deter- 
minant of the linear canonical transformation is 1 and that 9 and ip are two angles 
varying in an interval of length 2ir. Since m and n are relatively prime, the Bezout 
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theorem states that this equation has a solution (uo,vo). There are an infinite num- 
ber of solutions which can be written (uq + kn,Va + km) with k £ Z. We denote by 
(u, v) one of these solutions. A choice of a couple (u, v) is associated to a choice of 
a particular basis of the homology group. 

The generating function F 2 of type 2 [ Arn89 associated to the canonical trans- 
formation is given by 

(70) F 2 = (uipi + vtp 2 )J + (ntpi + mip 2 )l , 

where J and X are the momenta conjugated respectively to 9 and ip. From the 
definition of F 2 , one deduces that 

J = mil — nl 2 
2 = ul 2 — vli ' 



(71) 



and that as expected J = J. We will drop the tilde in the rest of the proof. 
The angle 9 can therefore be written as 

(72) 9 = u arg(pi +iq 1 )+v arg(p 2 + i®) ■ 

Differentiating Eq. (|71ll with respect to time and using the Hamilton equations, 
one obtains that 

A u ( n ^ + glPi + flgi) ^ <rnn 3 + ^p 2 + §§ q2 ) 

(73) e= WTpJ + W+ti - 

which leads to 
(74) 

. u (nmh-nmR + m^ Pl +m^ qi ) v(nmh - nmR + n§£p 2 + njgg 2 ) 



j + TTl TTl - j 

The last step consists in rewriting this integral as an integral in the reduced phase 
space Pj. The independent part of gives a continuous contribution on the line 
C denoted Oq. Since tti — 2mnn 2 , one finally obtains 

hu hv diri 



(75) &(hj)= + :]— + &o(h,j). 

Jtv- J + ""I ""1 J ^2 

For r the proof is straightforward and similar to the one of lemma [3] □ 
Lemma 7. The complex discriminant locus A near the origin is given by 

(76) A = I h = ± / (-l)-2" + "m"n"j'" + " 



(m+n) m + T) 

In the limit h —* 0, j < fixed, the polynomial Q as a function of x has n roots Xk 
whose leading term is 

(77) Xk = -, : — i—e « 

v > (_2j)m/n 

with k £ {0, 1, ■ • • ,n— 1}, the other roots having a finite limit different from zero. 

In the limit h — ► 0, j > fixed, the polynomial Q as a function of x' has m roots 
x' k whose leading term is 

u2/m 

(78) x' k = . . . e 
v ' k (2j)"/ m 

with k £ {0, 1, • • • , m — 1}, the other roots having a finite limit different from zero. 



1 • r i. 
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Remark 5. We notice that the case j > and j < give the same expressions but 
with m and n interchanged. 

Proof Let us assume that j < 0. As for the 1 : — n resonance, we calculate the 
complex discriminant locus of Q near the origin as a function of x. The construction 
of the Newton polyhedron gives the principal part Qn of Q 

(79) Q N (x) = (x-2j) m x n -h 2 . 

Simple algebra leads to the discriminant locus A. 

In the limit h — > 0, j < fixed, the construction of the Newton polygon of Q 
leads to the following equation for the roots of the principal part of Q 

(80) (-2j) m x n = h 2 , 

and to the n roots Xk of Eq. (|77[) . Exchanging the role of m and n and taking 
j > 0, a similar proof gives the roots of Eq. (|78j) .D 

Having established lemmas El [6] and [7] required, we can pass to the proof of 
proposition [3l 

Proof Since the line C is crossed two times by the loop T, the monodromy matrix 
M has two contributions denoted M_ for j < and M + for j > 0. We first consider 
the case j < 0. The case j > will be deduced from the calculation for j < by 
exchanging the role of m and n. The calculation is based on the analysis of the 
roots of the polynomial Q. Q has iV roots denoted Xi with i 6 {0, 1, • • • N — 1}. 
The roots x n , ■ ■ ■ xn-i are of order 1. A simple analysis of the polynomial Q leads 
to the following relations 



(81) ^ 



nf=o 1 ^ = ^ 2 (-i) JV+1 

1 li=0 — (-2j) m V -V 



Following the same steps as in the proof for the resonance 1 : — n, one arrives to 

p * dx 



(82) Q(h,j)-Q (h,j)~uh 



where ( J^)l"/„ < k < 1. For h > 0, Eq. J82]) simplifies into 

(83) e(/»,i)-e„(M~ 



+ OC 



x\/x n — 1 
We finally obtain that 

(84) lim = 6(/i, j) - o (/i, j) = ±u- . 

h— i-o±,j<o n 

A similar proof leads to 

(85) lim =Q(h,j)-Q (h,j)=±v- . 

h^>o±,j>o m 

The calculation of r uses the same arguments and shows that r is continuous on 
C. 

One then deduces that the matrices M_ and M + are respectively given by 

I 

_H 1 



(86) M_ = 
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and 

(87) M + = ( I ° ) , 

where we have taken into account for M + the fact that the line C is crossed from 
h < to h > 0. The total monodromy matrix is given by the product of the 
matrices M_ and Af + 

(88) m=m ;)(» ;)=(_! ;) , 

where the relation mu — nv = 1 has been used.D 

3. Extension to the complex domain 

3.1. Idea of the method. In this section, we reformulate the notion of fractional 
hamiltonian monodromy by deforming the loop T close to the line C, such that 
it bypasses the line C through the complex domain. The starting point of the 
complex approach is given by the expressions of the functions and r as real one- 
dimensional integrals (see for instance lemma [TJ Eqs. (fit))) and for the l:-2 
resonance case). If we consider the complexified variables h, j and 7Ti then Q(h,j) 
and r(h,j) can be interpreted as complex integrals on a line of the complex plane 
7Ti. Furthermore, these two functions can be rewritten as integrals of rational 1- 
forms over a cycle S on the Riemann surface defined by 7r| = Q(ni). This allows 
the use of topological properties of the Riemann surfaces. 

More precisely, we first introduce a Riemann surface constructed from the energy- 
momentum map F and we determine the Gauss-Manin connection of a complex 
semi-circle Tc bypassing the line C of singularities. As explained below, we can 
then construct the extension to the complex domain of the functions and r along 
Tq- The monodromy matrix is determined as in the real approach by the variation 
of the functions 8 and t along a loop T around the origin. The regularizations 
of the functions and r to cross the line C are replaced by their continuations 
along Tc- We obtain the fractional monodromy matrix by letting the radius of the 
complex semi-circle Tc tend to 0. 

3.2. Extension to the complex domain of the l:-2 resonant system. 

Remark 6. All that has been established in Sec. \2.1\ for the real approach can 
be done exactly in the same way for the complexified phase space T*C 2 where 
( J, tti, 7T2, ttz) € C 4 , except for the fact that there is no restriction on the values 
of J and 7Ti. The quotient is taken to be C* ~ S 1 x R*. The reduction in the 
complex approach is described in appendix^ The manifold C / S 1 x K* has real 
dimension 6. 

We begin by recalling some of the basic elements of the theory of complex alge- 
braic curves which will be used throughout this section (see |Kir93j for a compre- 
hensive introduction). 

Let P : C 2 — > C be a polynomial function. The fibers P^ 1 (z) are complex 
algebraic curves (hence two-dimensional real surfaces). There exists a finite set 
£ C C (essentially critical values of P) such that all fibers P _1 (z), z G C \ S, 
look alike. Moreover, the mapping P : C 2 \ (P^ 1 (E)) — > C \ S is a locally trivial 
fibration. Hence, given any path 7 : [xoj^i] — > C \ S, one can identify the fibers 
P -1 (7(2:0)) and P _1 (7(zi)). This identification is not unique but induces a unique 
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h=-h h=-ih„ h=h„ 

o oo 

Figure 6. Transport of the cycle 6 along a loop around the line 
C. The loop is an arc of a circle of radius ho > 0. The large solid 
lines represent arbitrary branch cuts of the Riemann surfaces and 
the full dots ramification points (see text). The parts in solid and 
dashed lines of the loop respectively lie in the upper and lower 
leaves of the Riemann surface. 

identification between the homology groups of the two fibers. For z$ £ C \ S, a 
cycle 8(zo) £ i?i(P fa)) can be transported along any path 7 in C \ S starting 
at zq giving thus a family of cycles S(z). The transport depends only on the 
homotopy class of the path 7 in C\ S. This is the classical Gauss-Manin connection 
AGZV88, Z0IO6 . The Gauss-Manin monodromy is defined from the Gauss-Manin 
connection as an automorphism of the homology group Hi for a loop in C \ E. 

For the l:-2 resonance (see Sec. 12.21) . the energy-momentum map of Eqs. ((T|) 
can be treated in the formalism of complex algebraic curves even if the situation 
is more complicated than described above. The difficulty here lies in the fact that 
the fibration is defined only implicitly by 

(89) F c := irl - [fa - j)fa + jf - (h - efa - j 2 )f] = , 

i.e. tt| = Qfa, j, h). 

For fixed generic values of (h, j), Eq. (|89|) defines a torus from which two points at 
infinity have been deleted. The corresponding complex algebraic curve is schemat- 
ically represented in Fig. [6l This representation can be understood by solving 
Eq. with respect to 112- Generically, there are four values of 7Ti for which 

Qfa, j, h) = 0, giving each a single solution for 7T2 = 0. These points are the ram- 
ification points of the complex algebraic curve. For all other points 7Ti, there are 
two solutions 7r 2 of Eq. (|8"9")) represented by the two leaves in Fig. [S] We denote by 
Ffj h s the set of points {tt 1i tt 2 ) £ C 2 such that fa, ir 2 ,j, h) verify Eq. l[8"9"]). Note 
that for two different values (h, j) ^ (h',j') the fibers F9- w and F9., h ,, intersect. 
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Nevertheless, it is possible to generalize the Gauss-Manin connection to this case. 
Consider for that the mapping G : C 4 — > C 3 , given by 

(90) G(7r 1 ,7r 2 ,/i,j) = (7r 2 2 -Q(7r 1 ,/i,j),/i,j) . 

The mapping G defines a fibration on the complement C 3 \S of the set E £ C 3 where 
its rank is not maximal. We use here the Ehresmann fibration theorem [Wolf64J. 
We take as basis of our fibration denoted BeC 2 the set ({0} x C 2 ) C (C 3 \ S) which 
is viewed as a set in the (h, j)-space C 2 . As in the real approach, we introduce the 
variable x = j +tti and we set y = Tt2- The singular locus of the fibration is the 
set of points (h,j) £ C 2 where the polynomial Q has multiple complex roots. This 
set is given by Eq. (|T5|) and is displayed for j 6 1 in Fig. [8] We notice that for 
(h, j) = (0, 0) the fibration is singular and that the singularity (x = 0, y = 0) is not 
of Morse type. 

In the reduced phase space Pj, the original real torus projects to a cycle 5(h,j) 
delimited by and 7T+ (see for instance Fig. HJ. 71 T an d "i are the two largest 
real roots of the polynomial Q as a function of ir± . Returning back to the Riemann 
surface and following notations of Eqs. (f!8|) . the roots of Q, which are simple 
ramification points of the Riemann surface, are denoted Xk (fc = 1, • • • 4). One can 
introduce cuts along the segment £2X3 an d along a simple curve joining x\ and X4 
and avoiding the segment x 2 x 3 . For (h,j) € lZ reg , the cycle S(h,j) is represented 
by the real oval between the two largest real ramification points which correspond 
respectively to J+Trf and j + Tit . To be coherent with the real approach, this cycle 
is oriented from x = j + to x = j + in the upper leaf and from x — j + 
to x = j + Hi in the lower one. All these notations are displayed in Fig. [6] 

3.3. Computation of fractional monodromy from the Gauss-Manin mon- 
odromy. We pursue in this section the construction for the l:-2 resonant system 
to arrive to the computation of fractional monodromy at the end of the section. Let 
r be a loop around the origin. We recall that the computation of the monodromy 
matrix associated to T is based on the difference of the values of the functions 
and r at each side of C as h goes to 0. The goal here is to compute the variations 
of these functions using their extensions to the complex domain near C. 

Definition 2. Starting with the result of lemma\]\ we introduce the complex con- 
tinuation of the functions and r defined by 



(91) 



■j) xy 
T {h,j) = Tii Is(h,j) 



where (h, j) 6 B and y 2 = (x — x\){x — x 2 ){x — x%){x — X4). The positive and 
negative determinations of the square root y have been chosen respectively for the 
upper and the lower leaves of the Riemann surface. We have added a factor i in 
the definition of and r to coincide with the real case. 

We locally deform in a neighborhood of the line C the loop T. We denote by 
r<7 this complex deformation and by Tr the rest of the loop. T, Tc and Tr are 
represented in Figs. [8] and [9] 

Remark 7. In this work, we have considered deformations of the real loop in the 
half-plane 3 [ft-] > but they could be equivalently done in the half-plane 3 [h] < 0. 
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Figure 7. Definition of the cycle So and Si for the l:-2 resonance. 
The position of the pole of 8 is represented by a cross. 

The bypass To (ho) is a semi-circle of radius ho in a plane with jo < fixed 
around the line C. The corresponding real path completing Tc(ho) is denoted 



^R(h-o)- From Eq. (| 1 5|) and for a small jo, one deduces that if ho < J -^rJo then 



X2 and X3 exchange their positions along Tc{ho) whereas if ho > \J ~^3o then 
3 ramification points of the Riemann surface (x\, X2 and X3) move and exchange 
their positions. The change of the ramification points is displayed in Fig. 1101 It 
can also be deduced from the asymptotic expansions of the roots of Q [see Eqs. 
(|18p]. This can be seen by parameterizing the loop Tc(ho) as 



where t G [0,7r]. There is thus a qualitative difference of the result depending on 
the value of ho- ho must be chosen sufficiently small to be in the first case since for 
each fixed jo we are interested in the limit h — * 0. 

The family of cycles S(h,j) can be obtained by transport of the cycle S(h ,jo) 
along r^. Examining Fig. [Till one sees that the lines of singularities crossed by 
Tr have no incidence on the ramification points defining the cycle S. The parallel 
transport of 5 along Tc(ho) is given by the change of the ramification points along 
T c (h Q ). From the Picard-Lefshetz theory |Zol06| IAGZV88j (see Figs. 0, one 
can show that the cycle S(ho,jo) when transported along the bypass Tc(ho) is 
transformed into S(—ho,jo) + So(— ho, jo) where So is a vanishing cycle around the 
ramification points X2 and X3 of the fiber F^_ ho . X2 and X3 are defined by Eqs. 
(Unj) - For the position of the cuts of Figs. the cycle S {—h ,jo) is composed of 
a path from x% to X2 on the upper leaf followed by the lift of the same path to 
the lower leaf run in the opposite direction. The vanishing cycles do and <5i are 
represented in Fig. The cycle Si will be used in lemma [8] Note the different 
choices of cuts for the Riemann surface between Figs. [6] and Fig. [7l Hence, an 
abelian integral I{ho,jo) — Js(h j ) 10 becomes after going once around the semi- 
circle Tc(ho) the sum I(—h ,j ) + fs (-h ja) 10 wn i cn gi yes the variation of the 
function / over Tc(ho). This latter remark can be applied to the functions 6 and 
r. However, due to the presence of a pole at x — for the function <d, the cycle 
So has to be positioned with respect to x = 0. The counterclockwise turning of the 
ramification points X2 and X3 implies that the cycle $o avoids the singularity x = 
from above. 




(92) 
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Figure 8. Complex discriminant locus A (solid lines) of the 
energy-momentum map of Eq. fl} for j £ 1 and h S C. The grey 
plane corresponds to the real bifurcation diagram. The dashed 
lines represent the loop T locally deformed near C to the complex 
domain. The arc of circle is in a complex h — plane with j fixed. 




Figure 9. Decomposition of the loop T into the paths Tc and Tr. 

From the analysis of the behavior of the cycle S along T, we can determine the 
variations of the functions 8 and r along this loop. These variations are denoted 
AOr(fto) and Arp(ho) and defined as follows 

Definition 3. For the energy-momentum map of Eq. Q), A0r(/io) o,nd Arr(/io) 
are given in the complex approach by 

AQ r (h ) = A&r G (h ) + AGr R (h ) 
Ar r (/i ) = ATr c {h ) + Arr R (/i ) 



(93) 



A simple calculation allows to simplify the expressions of AOr(/io) and Att(/io)- 
Since 

,~as /. ho f dx —ho f dx 

(94) AQ R (h ) = — - - — / - , 

Zl£ Js(h a ,j ) X V l%e Js(-h ,j ) X V 

and 

ho f dx —ho f dx ho f dx 

6(-h ,j ) xy 2 «£ J5 (-h ,j ) xy 2ie J S (h ,j ) xy 
one deduces for j < fixed that 

A8 r (^o) = -P- L , h ■ n — 

&Mho) = m I So{ - ho , jo) 



(95) Ae c ( ho ) = ^[ d JL + =^f 

2«e Js(-h ,i ) xy 2ie J So 



(96) 
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Figure 10. Schematic representation of three of the roots of the 
polynomial Q as a function of h and j. These roots are the roots of 
the principal part Qn of Q defined by Eq. (fTi|) . The polynomial Q 
has another root larger in module which is not represented here as 
it undergoes no bifurcation. The small inserts depict the graph of 
Qn as a function of x for different values of h and j. The position 
of the inserts gives the corresponding values of h and j. The solid 
and dashed lines are lines of singularities of A. The dashed lines 
do not belong to the real bifurcation diagram (see Fig. [5]). 



We introduce for the function the following quantities 

r Ae r = iim M Ae r (M 

(97) I A6 rc -lim w Ae rc (^o) , 
[ A6 r „ - lim^^o A8 rn (^o) 

and the same for r. 

Different results have to be established before computing the monodromy matrix. 
The variation of O around C can be viewed as a residue. 

Lemma 8. The sum of the integrals of the 1-form over So and S\ is inde- 

pendent of h Q and j and equal to 

(98) !*f d -l + ^L( ^ = 2n. 

2lE JSo(ho,jo) X V 2lE Js^hodo) X V 

In this equation, ho is taken sufficiently small and jo < 0. 

Proof We use the notations of Fig. [7] The union of $o and <5i corresponds to two 
loops around the pole x = lying respectively in the upper and the lower leaves 
of the Riemann surface. The orientation of these two loops is on the lower leaf the 
opposite to the one on the upper leaf. The same applies to the determination of 
the square root y for the two leaves of the Riemann surface. Hence, the sum of 
the left hand-side of Eq. (1081) is given by two times the residue of the 1-form — at 

xy 

x = 0. One deduces that 

ho f dx ho f dx ho , 1 

99 -f- / — + 7T- — = -^vriRes — ,a; = . 

2ie Js (ho,jo) x y 2i£ Ai(/«,jo) x v 2i£ x y 
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Simple algebra leads to 

1 £ 

(100) Res( — ,x = 0) = — , 

xy h 

which completes the proof. □ 

Lemma 9. The variations of the functions and r along T are 

(101) / A9 r ^ 

[ W > { Ar r - ' 

Proof We use Eqs. (|96| and we calculate these two quantities in the limit ho — > 
and jo < fixed. The asymptotic expansion of the roots x% for h — > and j < 
fixed is given by Eqs. (fTO]) [see Ref. |ECS07j for the explicit computation]. 

We begin by the computation of AOr- We first notice that x-i and x% are real. 
do can thus be viewed as a real loop which is locally deformed in a neighborhood 
of x = to avoid the pole in x = 0. Since the path So is oriented in the opposite 
direction and the sign of y is the opposite in the lower leaf with respect to the upper 
leaf, it is straightforward to see that the contributions of the upper and lower leaves 
of the Riemann surface coincide. AOr can be written as follows 

r -h f X2 dx h 1 , 1 o 

(102) A8 r = lim \PV / — + — -27riRes( — , x = 0) , 

ho->o is J X3 xy ie 2 xy 

where PV denotes the principal value of the integral. The introduction of the 
principal value is due to the presence of the pole at x = 0. The residue is calculated 
with the positive determination of the square root y. The factor ^= in front of the 
residue corresponds to the fact that the integral is taken on a semi-circle which is 
oriented in a clockwise manner. The contribution of the residue term to AOr is 
equal to ir. Following Ref. [ECS07], the computation of the principal value term 
can be done by using elliptic integrals. It can be shown that this term is zero. 

Using the same arguments, we can deduce that Arr = since r has no singularity 
along the real segment X2X3. □ 

In the limit ho — > 0, So and 81 play a symmetrical role for O. More precisely, we 
have 

Corollary 1. The integrals of the 1-form j^^j over So and S\ are given by 

(103) lim^/ ^=l im J^/ ^ =7r . 

h B ^o 2ie J 5o (h j ) ^0-0 2ie J 5l {h Q ,j ) xy 

Proof The proof of lemma [5] has already shown that 

ho f dx 
(104 lim — / — = 7T . 

h ^o 2ie J So (h ,j ) xy 

We conclude for 8\ by using lemma [HID 

Remark 8. We remark that the complex continuations of the functions O and r 
along Tc have replaced the regularizations of these functions in the real approach. 
The semi-circle Tc is taken to be asymptotic in order for AOr and Arr to be inde- 
pendent of h and j and to recover the topological character of fractional monodromy. 
In contrast, if we consider a loop around the line C then the variation of O along 
this loop is topological as it is calculated from a residue (see lemma[8ty. 
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From lemma [9j we can finally conclude by the following proposition. 
Proposition 4. The monodromy matrix M for the loop T is given by 



(105) 




Proof We use lemma [9] and the fact that the monodromy matrix is determined 
by the variations A0r and Arr-D 

Remark 9. In the computation of lemma\^ we have taken arbitrary j, which shows 
the topological character of the definition of fractional monodromy. 

We note that the real and the complex approach can be related by the following 
corollary. 

Corollary 2. A6 C = and hence AQ R = A6 r 



Proof We apply the Jordan's lemma to show that AOc = 0.D 

The monodromy matrix in the real approach (resp. complex approach) is given 
by A6_r and At/j (resp. AOr and Att). The corollary [2] shows that these two 
approaches are equivalent. 

3.4. Generalization to 1 : — n resonance. All the arguments used for construct- 
ing the extension to the complex domain of l:-2 resonant systems can be generalized 
to 1 : — n and m : —n resonant systems. As for the real approach, we consider the 
family T of energy- momentum maps introduced in Eqs. (I35| . 

The Riemann surface is defined from the relation deduced from Eqs. (|34[) and 

(25) 

(106) y 2 = x n (x - 2j) -[h- R(x — j,j)} 2 . 

The discriminant locus A of these systems is given locally by Eq. ([ST]) . Note that 
this locus is qualitatively different according to the parity of n. More precisely, the 
real lines of singularities for n odd become purely imaginary for n even. This does 
not change the discussion of this section. Following the preceding case, we consider 
a loop r around the origin which decomposes into a complex semi-circle Tq around 
the line C and a real path Tr. Along Yc, one sees by using the expansion of the 
roots of Q (lemma |4} that for j < 0, n roots in the variable x turn asymptotically 
around the origin by an angle The other roots stay fixed to first order in h. 
Fig. fTTI illustrates this point for the resonance l:-3. Note that only the principal 
part of the roots are exactly exchanged among each other. Following the change 
of the ramification points, we can transport the cycle S along To- S(h,j) is a real 
oval between xq and the ramification point of the Riemann surface associated to 
71"^ denoted x n . This cycle is oriented from xq to x n in the upper leaf and from x n 
to xq in the lower one. Figures Q1J] and [T2] display the transport of 5. Not all the 
ramification points are represented in Figs. 1131 and [141 The position of the cuts is 
arbitrary but indicates on which leaf of the Riemann surface a path lies. We can 
compare two Riemann surfaces if the cuts of the two surfaces are the same. Since 
for the l:-n resonance the cuts move after a semi-circle, one has to modify the cuts 
to recover the initial choice of cuts. An example of this deformation is given in 
Figs, m 
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9t(x) 



Figure 1 1 . Evolution of three of the roots (see text) of the poly- 
nomial Q from the resonance l:-3 and for a complex semi-circle 
going for ho > to — ho < (jo < fixed). Numerical values are 
taken to be ho — 0.005 and jo = — 1- The radius of the semi-circle 
is ho- The energy-momentum map is given by Eqs. (|36| . The full 
and open dots represent respectively the roots for the starting and 
the ending points of the path. The dashed line is a circle of radius 

h 2/3 

(_2j)i/3 wmcn corresponds to the leading term of the expansion of 
the roots as ho — > 0. 



Figure 12. Definition of the cycles So, 5% and 62 for the l:-3 
resonance. The cross indicates the position of the pole of 0. 



We see that after the semi-circle Tc(ho) around C, 5(ho,jo) is transformed into 
S(—ho, jo) + So(—ho, jo). So is composed of a path from xo to x n -i in the upper leaf 
and of the lift of the same path in the lower leaf but run in the opposite direction. 

Definition 4. We define the complex extension of the functions and r as 



5 



1 




(107) 
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Figure 13. Transport of the cycle 5 along a semi-circle around 
the line C . The radius of the semi-circle is Hq. The last three 
figures arc equivalent but with different cuts. The cross indicates 
the position of the pole of 0. 





h=h 



h=-h. 




h=h 



h=-h 



Figure 14. Same as Fig. [13]but for 3 semi-circles around the line C. 



where the positive determination of the square root y is associated to the upper leaf 
and 

1 f x n ~ x dx 
2n 



(108) 



>s(h,j) y 

Qq is the complex extension of the real function Qq introduced in Eq. ^By . 0o has 
a trivial contribution to the monodromy matrix. 

We define the variations AOr(h ) and Arr(/io) of the functions and r along T 
as in the definition [3] for the l:-2 resonance. These variations are given by integrals 
over the cycle So- A6r and Arr denote the limits of these variations as ho — > 0. 
They are given by 

A9 r =lim /lo ^ Ae r (/io) 
Ar r = hm,, ^o Ar r [ho) 

Equivalent lemmas to lemmas [5] and [S] for the 1 : — 2 resonance can be established 
for the 1 : — n resonance. We only state the final result. 



(109) 



Proposition 5. The monodromy matrix associated to the loop T is equal to 



(110) 



M = 



1 
-i 1 



Proof We compute the different integrals in the limit ho 
A9r(ft-o) can be written as 

-ho f dx 

2 Js (-ho,jo) X V ' 



with jo < fixed. 



(Ill) 



A6 r (M 
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where So is the cycle between xo and x n -i which is oriented from xq to x n -i in the 
upper leaf and inversely in the lower leaf. In the limit ho — > 0, one deduces that 

r xa dx 



(112) A9 r = lim h [ 

h o^0 J 



- 1 X V Uk=o( x ~ x k) 



where the integral is taken along an arc of a circle from x n —\ to xq of radius („ 2 j) 1 /" 

h 2/n 

Using the change of variables x = f _° 11/ „ e %x , we obtain 



(113) A6 r = 



27r/ " d x 



m=o^ x -^ kln ) 



It is straightforward to check that the integrand is not modified by the translation 
X' = X + ~ ■ One then deduces that 

(114) A9 r = lim --Res[ — ~ H ° =,x = 0] . 

x y Uk=o( x ~ x k) 

Simple algebra finally gives 

(115) A6 r = — . 

n 

Similar arguments show that Arr = 0. We finally construct the monodromy matrix 
from the variations of 8 and r along T. □ 

We finish this section by presenting a complementary computation of fractional 
hamiltonian monodromy in the complex approach. The idea is here to use direct 
asymptotic computations to determine AOr- 

Lemma 10. For j < fixed and \ h\ — ► 0, we have the following asymptotic behavior 

(116) 9(h,j)-Qo(h,j)~ ~arctan[^] , 



where ko is a function of h such that lim^i^g fe^ 



0. 



Proof We proceed as in the real approach by using the asymptotic expansions of 
the roots Xk of the polynomial Q (see lemma H]). Following computations of the 
proof of proposition [2l we obtain 

P x " dx 

(117) e(h,j)-e (h,j)~- 



h 2 /" / T1 i=N-l, 7 



which transforms into 

(118) e(h,j)-e (h,j)~h 



k dx 



xJ(-2j)(x»-X,) 



The real function k fulfills lim|^|^o ^^P/n — = 0. Introducing the variable u such 

i 2/n 

that x = ^_ 2 -y/ n u, one arrives to 

(119) e(h,j)-e (h,j) 1 "' 



ly/u n - 1 
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Figure 15. Riemann surface of the function arctan. The solid 
and dashed lines respectively lie in the upper and the lower leaves. 
The ramification points correspond to the complex numbers i and 



where we have used the fact that \fh 2 = h. One finally obtains that 

o fc(-2j) 1 /" 



(120) G(h,j)-Q (h,j)~ [- arctan [Vu^T]] x 
which leads to 

(121) e(h,j)-Q (h,j)~ -arctan [^] , 

where k = k n / 2 ^2j.D 

The behavior of the function near the line C is thus related to the complex 
function arctan. Using the fact that 

1 1 + iz 

(122) arctan z = — ln[ ] , 

2i 1 — iz 

for z G C, we can construct the Riemann surface of this function. This surface has 
two leaves and a cut between the points z = i and z = —i. Figure [TBI displays this 
surface in the variable z. We recall that this function has a jump of it along a loop 
crossing once the cut. Using Eqs. (|57|) and ([58]) . one sees that the limits of O when 
h —* correspond respectively to the points i?+ and i?_ of the Riemann surface 
(see Fig. [15]). Since these two limit points are the initial and final points of Tr, 
we can associate to Tr the real path from R_ to R + of the Riemann surface of 
the function arctan. The complex continuation of the function along the small 
semi-circle Tq of the bifurcation diagram is associated to the big circle at infinity 
of Fig. [15] From Fig. Qjj] we thus recover that after a real loop T locally deformed 
to the complex domain to bypass the line C the function has a jump of — . This 
result is also coherent with the analysis performed with the variation of the cycles 
5. We denote by Sk the cycle between the ramification points Xk and Xk+i- The 
cycle is oriented from Xk to Xk+i in the upper leaf and from Xk+i to Xk in the 
lower one. Following the proof of proposition [5] it can be shown that if we compute 
along one of these cycles then we obtain asymptotically the same result. More 
precisely, we have 

h f dx 2t: 



(123) 

2 Js k (h,j) xy 

for j < fixed and h — > CU. We consider now r semi-circles around C. The cycle S 
is transformed into 6 + So — Si + S 2 — ■ ■ ■ + (— l) r ~ 1 5 r -i. This is displayed for r = 3 
in Figs. HU As expected, we thus see that after an even (resp. odd) number of 
semi-circles, the function has no jump (resp. a jump of ^) which corresponds 
to the behavior of the complex arctan. 



:S2 



D. SUGNY^P. MAR.DESIC 2 ,M. PELLETIER 2 , A. JEBRANE 2 ,H. R. JAUSLIN 1 



3.5. Generalization to m : — n resonance. We now study the m : —n resonant 
system with m > 1 and m and n relatively prime. We associate to such a system a 
Riemann surface which can be constructed in the variables x = it\ +j or x' = tti —j 
from 

(124) y 2 = x n (x - 2j) m -[h- R(x - j,j)f , 
or 

(125) y 2 = (x' + 2j) n x' m -[h- R(x' + j,j)] 2 . 

One passes from one representation to the other by the relation x = x' — 2j. In 
particular, the two surfaces have the same ramification points translated by 2j. 
Depending on the line of singularities considered, one or the other surface will be 
used, i.e., the surface in x for j < and the surface in x' for j > 0. We next recall 
that the discriminant locus A of such a system is given by Eqs. (|76|) . From the 
expansion of the roots of the polynomial Q in h = (lemma [7]) , we deduce that for 
j < (resp. j > 0), n (resp. m) roots in x (resp. in x') exchange their positions 
along a loop around the line C. We locally deform T in a neighborhood of the line C 
and we decompose this loop into four loops , Tr 2 , Td and Tc 2 where and 
Tc 2 are respectively two semi-circles around the lines of singularities (h = 0, j < 0) 
and (h = 0,j > 0). and Tji 2 complete the loop T respectively for h < and 
h > 0. We define the variations of the functions O and t along T as follows 

(126) { A6r = AGrc i + A6r ^ + A9r "! + A9r "- 
1 ' I Ar r = Arr Cl + Arr C2 + Arr fll + Arr fl2 

where the different variations are determined asymptotically, i.e., when the radii of 
the semi-circles C\ and Ci go to 0. A simple calculation then show that the study 
of the to : —77, resonance can be reduced to the study of two cases for j > and 
j < similar to the 1 : — tt, resonance. This allows to compute the jump of the 
function along T and to deduce the monodromy matrix for the real loop T. We 
recover the result of proposition [3] 

4. Conclusion 

In this paper, we have investigated the notion of Fractional Hamiltonian Mon- 
odromy in the l:-2, 1 : —n and m : —n resonant systems. We have discussed our 
asymptotic method to calculate the monodromy matrix in the real approach and 
we have proposed a definition of fractional monodromy using a complex extension 
of the bifurcation diagram. From this definition, we have recovered the results of 
the real approach. At this point, the question which naturally arises is the gener- 
alization of this concept to other types of singularities. A part of the answer could 
be given by applying the complex approach to a new generalization of standard 
monodromy, the bidromy which has been recently introduced in Ref. |SZ07j . 

Appendix A. Geometric construction 

A.l. l:-2 resonance. As stated in the introduction, the geometric construction 
of fractional monodromy generalizes the construction of standard monodromy. We 
propose a schematic representation that is slightly different from the original one 
proposed in Ref. NSZ06 . This representation has the advantage to be easily 
generalizable to l:-n resonant systems. We use a standard representation of a torus 
i.e. a rectangle whose edges are identified according to the arrows. The crossing 
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Figure 16. Schematic representation of the continuous transport 
of the basic cycles (/?i,/?2) when the singular line C is crossed. 71 
and 72 are respectively represented in dashed and solid lines. 
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Figure 17. Schematic representation of the basis cycles for h < 0. 
Representatives of [fli] and [2/y are respectively represented in 
dashed and solid lines. 

of the line C is displayed in Figs. [16] and [TTj A curled torus is represented by 
two rectangles glued along a horizontal edge and with a particular identification 
of the vertical edges (see Figs. HH h = 0). This construction is not limited to a 
given energy-momentum map but can be applied to any energy-momentum map 
having a bifurcation diagram with such a line of singularities. Let /3i and /?2 be 
the two cycles associated to the flows of Xi and X2 defined in Eqs. ([6]). They are 
the representatives of the classes of homology [Pi] and [P2] which form a basis of 
Hi(T 2 (h,j),1) where (h,j) £ TZ reg . Only a sub group of Hi can be transported 
continuously across the line C. We assume that a basis for this subgroup is given 
by [Pi] and [2/32] • 2/3 2 corresponds to the cycle P2 covered twice. For h > 0, as 
representatives of and {2f3 2 ] we consider the two cycles 71 = (3i and 72. The 
cycle 72 is the union of two cycles /?2 with starting points belonging to the same orbit 
of the flow of Xi but separate by an angle 9 = n. To make a link with analytical 
calculations, 8 is taken to be ir/2 in this case. 71 can be easily transported across 
C and it remains unchanged. To transport continuously 72, the two cycles forming 
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Figure 18. Same as Fig. [16] but for the l:-3 resonance. 



72 have to be connected in one point in h = and then merged to form only one 
cycle covered once for h < 0. Representatives of the basis of Hi(F~ 1 (h,j),'Z) for 
h < are given in Fig. [T7] note that = —tt/2 since 8 has a discontinuity of 
size 7r on C. Comparison of Fig. [16] (h < 0) and Fig. [17] leads to the conclusion 
that after one loop, the cycle 72 becomes a representative of the equivalence class 
[2/3 2 — Pi] and that the corresponding monodromy matrix written formally in the 
basis ([Pi], [P 2 ]) is given by 



(127) 



M = 



1 

-1/2 1 



A. 2. l:-n resonance. The geometric construction of Sec. lA.ll can be straightfor- 
wardly generalized to other resonances. We consider the resonance l:-3 but other 
resonances can be treated along the same lines. The 3-curled torus is represented 
in Fig. [18] (h = 0) by three rectangles glued along a common edge. Note also the 
particular identification of the vertical edges. Following notations of Sec. IA.lt a 
basis of the subgroup of Hi(F~ 1 (h,j),'Z) which can be transported continuously 
across C is given by [Pi] and [3j8a], For h > 0, as a representative of [3/32] we 
choose three cycles P2 whose starting points belong to the same orbit of the flow 
of Xj with an angle 9 of 27r/3 between each other. We also assume that = tt/3 
in this case. After crossing C, this cycle belongs for h < to the equivalence class 
[3/32 — Pi] as shown by Fig. [TH] For h <0, 6 = — it/ 3 since the discontinuity of 6 
on C is equal to 27r/3. In the basis ([Pi], [P2]), the monodromy matrix M is given 

by 



(128) 



M = 



1 







-1/3 1 
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Appendix B. The semi-classical point of view 

We illustrate in this section the relation between classical monodromy and its 
semi-classical counterpart for 1 : — n and m : —n resonant systems. To our knowl- 
edge, this point has not been discussed up to now in the literature. We refer 
the reader to Refs. |Ngoc99[ INSZ06] for a rigorous definition of this semi-classical 
point of view. Here, we consider only the graphical representations of semi-classical 
monodromy as a pictorial illustration of classical monodromy. 

We consider two quantum differential operators J and H whose classical limits 
are the Hamiltonians J and H . These two operators commute i.e. [J, H] = 0. 
They thus have a system of common eigenfunctions belonging to L 2 (M. 2 , dqi A dqz). 
The corresponding eigenvalues form the quantum joint spectrum of the energy- 
momentum map F which is a 2-dimcnsional lattice of points. We can construct the 
quantum-classical bifurcation diagram of F by superimposing both the quantum 
joint spectrum and the classical bifurcation diagram. This has been done for the 
1 : — n and m : —n resonant systems in Figs. 1201 and[2"Tl Note that the quantum 
joint spectrum is defined for a given value H viewed here as a parameter. Using 
EBK quantification rules, we also introduce the semi-classical joint spectrum which 
is defined as the set of points (h,j) £ IZreg where the numbers m and given by 
the relations 

(129) h{m + ^) = j> pdq , 

are integers. In Eq. (|129[) . pdq is the Liouville 1-form and a* the Maslov index 
associated to the cycle where ([71], [72]) is a basis of Hi(T 2 (h,j),Z), We re- 
mark that the semi-classical lattice differs from the quantum lattice by o(K) which 
is irrelevant in our study. We also point out that the EBK quantification rules 
are not valid near the line of singularities C and have to be replaced by singular 
rules |dVN03j . In contrast, the quantum joint spectrum gives a smooth transition 
of the crossing C. Moreover, locally around a regular value of 1Z, this lattice is 
regular in the sense that there exists a map which sends this lattice to HZ 2 as h 
tends to zero. A systematic approach has been developed to check the regularity 
of the global spectrum. The method consists in taking a cell i.e. a quadrilateral 
whose vertices lie on the points of the lattice, transporting continuously this cell 
along a loop T and comparing the final cell with the initial one. If the two cells 
are different then the system has quantum monodromy. The rotation matrix which 
sends the initial cell to the final one is the quantum monodromy matrix Mq . More 
precisely, if the cell is supported by the two vectors (twijtua) and if after a loop 
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Figure 20. Semi-classical bifurcation diagram for the resonance 
1 : —3 and the energy-momentum map of Eqs. (|36p . The line of 
singularities C is represented by a solid line. The open dot indicates 
the position of the origin of the bifurcation diagram. The final cell 
after a counterclockwise closed loop around the origin is depicted 
in dashed lines. 



these vectors are transformed into (w' 1 ,w' 2 ) then Mq is defined by the relation 



(130) 



( < " 


1 =M Q 




V W 2 J 




V W 2 J 



Using the semi-classical joint spectrum, it can be shown that Mq = (ML) -1 
where Mqi is the classical monodromy matrix. We recall that when the system 
has fractional monodromy the size of the cell has to be increased to cross the line 
of singularities. In a way analogous to the classical simple cell cannot be 

transported continuously across the line C. The multiple cell then becomes the 
basic cell which is equivalent to consider only a sublattice of the original lattice. 
Here, the vertical lines of the quantum bifurcation diagram which are parallel to w\ 
are labeled by n\ , the quantum number associated to <p j which is a global quantum 
number. The cells are thus multiplied in the other direction. For instance, for the 
resonance l:-3 the size is multiplied by 3. This point is displayed in Fig. [20] For a 
resonance m : —n, due to the form of the monodromy matrix the size is increased 
from 1 to ran which explains why we have considered cells of size 6 in Fig. [5l] for 
the 2:-3 resonance. Examination of Figs. [20l and [2T1 shows that in the basis {w\ , w-x) 
the quantum monodromy matrices are respectively equal to 

<m, ( ; f ) , 

for the l:-3 resonance and 

(I-, ( \ r ) . 

for the 2:-3 resonance. The energy-momentum map used for the 2:-3 resonance is 



FRACTIONAL HAMILTONIAN MONODROMY FROM A GAUSS-MANIN MONODROMY 37 





r — ^— — . • * . • 
• . • 1 . " . 


• I 


1 — 1 




— I — t — '. — s — ._i 


• 1 

■ ■ * 


^""-^^ 1 





j 



FIGURE 21. Same as Fig. but for the 2:-3 resonance. The 
corresponding energy- momentum map is given by Eqs. p33|) . 



defined by 

(133) F 



J 

7T3 + (VTI + J) (TTI - J) 



Appendix C. Reduction in the complex approach 

We detail in this section the reduction in the complex approach for an m : — n 
resonant system. This reduction is different from the reduction used in the real 
approach which is associated to the S^-action of the flow of the Hamiltonian J. 
Both reductions leave invariant the polynomials (J, tti, 7T2, ^3) defined by Eqs. (|33[) . 

We start from the complexified phase space T*C 2 with (pi, qi,p2, 92) G C 4 . The 
reduction is based on two SO(2,C) actions $1 and $2. We recall that a matrix 
R G 5*0(2, C) is a 2 x 2 matrix which reads 



(134) 



a —b 
b a 



where (a, b) £ C 2 and a 2 + b 2 = 1. $ fe (k e {1,2}) is a map from 50(2,C 
C 2 which associates to each couple (7?, (qk,Pk)) the point of coordinates 



(135) 



Ofc -b k \( qk \ _ f a k qk - b k Pk 
bk a k J \ Pk J ~ \ b k qk + a k Pk 



We next introduce new coordinates which can be written as follows 



(136) 



m = 91 - wi 

£1 = qi + ipi 
m=Q2- ip2 
6 = 92 + iP2 



Under the action of these new coordinates transform into XkS,k and X7 r]k 
where Xk = ak + ibk- Finally, simple algebra shows that the complexified invariant 
polynomials ( J, tt±, 7r2, ir^) are invariant under the conjoint action of $1 and $2 if 
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A™ A™ = 1. One can conclude that the reduction is associated to a C*-action $ 
from C* x C 4 to C 4 defined as follows 

(137) $: (A,(6,r?i,6,%)) - (A ro £i, >T m m, A _n &, A"r? 2 ) . 
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